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Exact-exchange energy density and energy density of a semilocal density functional approximation 
are two key ingredients for modeling the static correlation, a strongly nonlocal functional of the 
electron density, through a local hybrid functional. Because energy densities are not uniquely 
defined, the conventional (Slater) exact-exchange energy density g"''^'^""^) jg jjot necessarily well- 
suited for local mixing with a given semilocal approximation. We show how to transform g™^'^""^) jjj 
order to make it compatible with an arbitrary semilocal density functional, taking the nonempirical 
meta-generalized gradient approximation of Tao, Perdew, Staroverov, and Scuseria (TPSS) as an 
example. Our additive gauge transformation function integrates to zero, satisfies exact constraints, 
and is most important where the density is dominated by a single orbital shape. We show that, 
as expected, the difference between semilocal and exact-exchange energy densities becomes more 
negative under bond stretching in HeJ and related systems. Our construction of g™''^""^' {-,y g, 
resolution-of-the-identity method requires uncontracted basis functions. 

PACS numbers: 31.15.Ew, 71. 15. Mb 



I. INTRODUCTION 



In Kohn-Sham density functional theory [l|, 0], the 
exchange-correlation (xc) energy i?xc must be approxi- 
mated as a functional of the electron spin-densities n| (r) 
and 7ij,(r). This functional can be always written as 



dv exc(r), 



(1) 



where exc(r) = n(r)exc(i") is the exchange-correlation 



energy density, n = rtt 



is the total electron den- 



sity, and Exc is the exchange-correlation energy per elec- 
tron. Approximations to excfr) can be constructed in 
a fairly systematic way 0, |j] by employing increas- 
ingly complex ingredients built from the Kohn-Sham or- 
bitals. Most of the existing exchange-correlation approx- 
imations use only ingredients found from the occupied 
Kohn-Sham orbitals at r or in an infinitesimal neighbor- 
hood of r, such as n„{v) — YlT'^ l0icr(r)P, Vn^[r), and 
To- = \ Y1T'' l^'/'i<T(i")Pi where a =t,i- Such function- 
als are called semilocal and include the local spin density 
approximation [lUl, Q , the generalized gradient approx- 
imation (GGA) [7f7and the meta-GGA ]!]■ 

Semilocal functionals are often accurate [1,0, [nnn, 

ITsj l but tend to make large errors for open systems of 
fluctuating electron number [3, 

EHElii, such as 



fragments connected by stretched bonds. This occurs 
because semilocal functionals respect the exact hole sum 
rule for a closed system but not for an open one of fluctu- 
ating electron number [iQ^ , where (after symmetry break- 
ing) the semilocal exchange typically overestimates the 
magnitude of the static correlation. We have argued [ll| 
that, in order to correct these errors, one needs to go 
beyond the semilocal approximation and incorporate a 
fully nonlocal ingredient, the exact-exchange (ex) energy 
density ex''(r) conventionally (conv) defined as 



^cx(conv) 



(r) 



E 

4,i 



(2) 



where 7o-(r, r') is the one-electron cr-spin density matrix 
of the Kohn-Sham reference system 



occ. 



0.,(r)0*o(r'). 



(3) 



In the Jacob's ladder classification of density functional 
approximations 0], functionals that employ e°''(r) are 
called hyper-GGAs. We have also argued that a hyper- 
GGA can simultaneously achieve good accuracy and sat- 
isfy important exact constraints if the exact-exchange en- 
ergy density is combined with a semilocal (si) exchange- 
correlation in a so-called local hybrid (Ih) functional 



-|-[l-a(r)](e^'-er) + e^ 



(4) 



2 



where < a(r) < 1 is the position-dependent mixing 
function. If a(r) — const, Eq. (g]) reduces to a global 
hybrid (gh) functional [l^l • The general local hybrid form 
was suggested by Cruz et al. [2l| as early as 1998, but 
specific forms of a(r) were not proposed until later [3l.l2^. 
The fundamental physical justification for local hybrids 
has been advanced only recently [l^. The local hybrid 
approach is, of course, not the only way of attacking the 
static correlation problem. Other distinct approaches are 
being actively pursued [H, . 

When a(r) of Eq. ^ tends to 1 in the high-density 
limit, the local hybrid functional uses full exact ex- 
change and treats correlation as the sum of two parts, 
the static (long-range, left-right) and dynamic (short- 
range) correlation. The dynamic correlation is relatively 
easy to model by a semilocal correlation functional e^' (r) . 
The static correlation is represented by the difference 
[e^'(r) — e°'*^(r)] weighted by a position-dependent func- 
tion [1 — a(r)]. This form is motivated by evidence that 
some (typically more than 100%) of the static correla- 
tion is already contained in semilocal exchange approxi- 
mations 19, 25, 26, 27], but not in e^''(r). 

Any proposal for a practical local hybrid functional 
must deal with the fact that, while the total energy is 
measurable, physical, and unique, the energy density is 
not. As a result, an arbitrary function G{r) that has a 
dimension of energy per volume and integrates to zero 
can be added to any energy density of a global hybrid 
functional with no effect on the total energy. In contrast, 
addition of G(r) to e™(r) or e^{r) in a local hybrid of 
Eq. (Ill will affect the total energy because the exact- 
exchange energy density here is weighted locally. 

While there is no "most correct" choice for the xc- 
energy density, there is indeed a conventional choice, 
which for exchange is Eq. However, the conventional 
exact exchange energy density does not have a second- 
order gradient expansion [28] and so is not the most natu- 
ral choice for density-functional approximation. Also, the 
exchange hole associated with g^''('=°'™) jg highly delocal- 
ized, which makes it very difficult to model with semilocal 
functionals. Standard functionals are at most designed 
to recover the conventional (or any other) exchange en- 
ergy density to zeroth-order in the density gradients, i.e., 
for uniform electron densities only. 

In the local hybrids proposed to date [13, H^, [13, [Ml, 
[3^ . (r) is taken as the integrand of the semilocal func- 
tional as written, while e"'' is taken as g™('=°'^'^) ^ This 
choice is not necessarily the one best suited for modeling 
the static correlation by the difference (e^ — e"'^). More- 
over, the very idea of attaching physical significance to 
the difference (e^ — e^^^) requires that both e^' and e"'^ be 
defined with respect to some common reference or gauge. 

The choice of the gauge itself is a matter of convention. 
One such choice is based on the Levy-Perdew virial rela- 
tion ^3] ■ Burke et al. have pointed out that virial ex- 
change energy densities e^"{r) = — n(r)r ■ Vz;x(r), where 
Wx(r) — dE^/5n{r), are unique for any given functional. 
However, the virial energy density depends on the choice 



of origin of r and has other undesirable properties. Fur- 
thermore, for the exact-exchange energy density, this ap- 
proach requires constructing the optimized effective po- 
tential [Sg, [13] (OEP), a procedure that is problematic 
in finite basis sets (371 [ssl [39j . Burke et al. have also 
proposed [U and investigated [13] the "unambiguous" 
exchange- (correlation) energy density which is uniquely 
determined by the corresponding energy functional via 
the exchange- (correlation) potential and the Helmholtz 
theorem. This "unambiguous" exact-exchange energy 
density has all the desired properties but, like the virial 
energy density, requires construction of the OEP and, 
hence, is not very practical at present. 

In this work, we propose and implement two new, de- 
pendable methods in which e^{r) serves as the reference 
and e°'^(r) is "tuned" to the gauge of e^{r). The first 
of them, summarized in section III Al below, is the one 
we will use in a still-unpublished hyper-GGA [4l| based 
upon the ideas of Ref . [13] ■ 

II. THEORY 

For a slowly- varying electron density, the conventional 
exact-exchange energy density will be well approximated 
by local or semilocal density functionals, although (un- 
like the integrated exchange energy) it has no analytic 
gradient expansion (28l . |43 ]. Our idea for making 
compatible with a given is based on the observation 
that although the static correlation is generally quite 
large (comparable in magnitude to exchange), it is neg- 
ligible in compact closed systems, such as atoms with 
nondegenerate electron configurations. Therefore, 
should be close to e^' at each r in such systems. This 
is consistent with the fact that the conventional exact- 
exchange [13, [3, [13] or exchange-correlation [13] energy 
densities in compact closed systems can often be modeled 
very accurately using only semilocal ingredients. 

Hence, we will make e^'^ as close as possible to e^' in 
those systems where the static correlation is known to 
be small. This can be achieved by various means: for 
example, by adding to g°''('^°"^) g, term that integrates 
to zero. We say that the resulting exact-exchange en- 
ergy density is in the gauge of that particular semilocal 
exchange approximation and denote it by ex''^^'^ ■ To illus- 
trate this method, we will construct the exact-exchange 
energy density in the gauge of the meta-GGA of Tao, 
Perdew, Staroverov, and Scuseria (TPSS) 

A. Construction from the divergence of a vector 
field 

For use as ex'^(r) in Eq. (|3|), we construct the exact- 
exchange energy density in the gauge of a semilocal func- 
tional as follows. First we write 

gcx(sl) (j.) ^ gcx(conv) (j.) +Gir), (5) 
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where g^^('=°'^'^) (^j.) jg ^j^e conventional exact-exchange en- 
ergy density given by Eq. ^ and G(r) is the gauge trans- 
formation term to be determined, such that 



J drG(r) = 0. 



(6) 



Obviously, Eq. ^ leaves much freedom in choosing the 
analytic form of function G(r). The range of possibil- 
ities can be narrowed down by several physical consid- 
erations: a) e^r'^^'^r) should reproduce e^{r) in atoms 
as closely as possible; b) for use in a hyper-GGA, G(r) 
should contain only the hyper-GGA ingredients, i.e., 
"<t(i"), To- = sEr''' l^<^*'^(i")P' eTJ-''°""'\r), and, possi- 
bly, their derivatives; c) Cx'^*'''''' (r) should satisfy as many 
exact constraints as possible. 

We start the construction of G(r) for spin-unpolarized 
systems by noting that the integral of the divergence of 
any well-behaved rapidly decaying vector field F(r) is 
zero, that is, /drV • F(r) = 0. So we will take G(r) = 
V • F(r). The vector field F(r) itself will be chosen from 

the requirement that ex''^'^°"^'' (r) + G(r) satisfy the most 
basic properties of the exchange energy density: correct 
coordinate scaling, finiteness at the nucleus, etc. 

One particular form that meets these requirements is: 



G(r) = aV • 



n/e 



^2 



l + c(n/£3)2 V T 



-W 



Ve 



(7) 



where e(r) = — ^"'^^ (r), = |Vnp/8n is the von 
Weizsacker [l^l kinetic energy density for real orbitals, 
T = -\-Ti is the Kohn-Sham kinetic energy density, and 
a, 5, and c (c > 0) are adjustable parameters. Note that 
< /t < 1 [48]. 

The function G(r) of Eq. ([7]) has the following ex- 
act properties of the exact-exchange energy density in 
the conventional gauge (or coordinate-transformed as de- 
scribed in Sec. IIIB[) : 

(i) Correct uniform coordinate scaling. Under this 
transformation, the conventional exact-exchange energy 

4gCx(conv)^^^^ or, 

A^. The ingredients of G(r) be- 



density behaves [33| like g^'^^'^""^^ (j-) 



in shorthand, 



have like n ^ X 



3 ? 



A, 



A^ 



A^ V - A, so 



G\{r) = A'*G(Ar), which is the correct behavior. 

(ii) Correct nonuniform coordinate scaling 49] . Under 
this scaling, the density behaves like ri^(r) = Ar7,(Aa;, y, z) 
or, in shorthand, n ~ A. The other ingredients scale 
in the A ^ cx) limit like e A°, ^ A^, r - A^, 
V ~ A, so in this limit G^{x, y, z) = XG{Xx, y, z), which 
is the correct nonuniform coordinate scaling property of 
the exchange energy density. 

(iii) G(r) is finite everywhere. This is because e has no 
cusp at the nucleus 50], which ensures that V^e is finite. 
All other ingredients of G are also finite. 

(iv) G(r) vanishes for a uniform electron gas and, more 
generally, satisfies Eq. 

We also note that, at large r, the density decays expo- 
nentially, n ~ e"""", where a is a constant, /t 1, 



e 1/r, so the large-r behavior is G(r) ~ —dn/dr n, 
which is comparable to the —n/2r decay of g°^('^°'^^)_ 

The values of a, 6, and c are determined by fitting 
gcx(conv)^^^ ^ Q^^^ ^,1^^^^ ^j^g^g sl^TPSS. In doing 

so, we note that for one- and closed-shell two-electron 
(iso-orbital) densities r^/r = 1, so G(r) is fixed by the 
parameters a and c alone. We use two model-atom iso- 
orbital densities: the exact two-electron exponential den- 
sity n{r) — (2/7r)e~^'' and the two-electron cuspless den- 
sity n(r) = (l/27r)(l + 2r)e-'^\ In the case of sl=TPSS, 
the fit gives a = 0.015 and c = 0.04. The value 5 = 4 
is chosen to be an integer that gives the best fit to the 
TPSS exchange energy density for the 8-electron jcUium 
cluster with rs = 4 bohr. This choice ensures that G(r) 
is very small (of the order of V^°) for a slowly varying 
density, as it should be. Gauge corrections for semilocal 
functionals other than TPSS can be constructed simi- 
larly by assuming the same analytic form for G(r) and 
refitting the parameters a, h, and c. 

While we cannot rule out that there exists a simpler 
function G(r) that satisfies exact constraints (i)~(iv), we 
can point out that many obvious candidates definitely fail 
to do so. For example, the function V^n^/'^, motivated 
by the work of Cancio and Chou [46| , correctly integrates 
to zero and has the correct uniform scaling property, but 
diverges at the nucleus and does not have the proper 
nonuniform scaling property. 

For a partly or fully spin-polarized system, the gauge 
correction becomes the sum of same-spin contributions 
= Y^a^<^i'^)- deduce the form of Go-(r) we use 
the spin scaling relation [sij : 

E.[n^,n^] = ^E^[2n^] + ^E^[2ni], (8) 

which also holds for exchange energy densities. Applying 
Eq. dH]) to G(r), we write 



G(K,nJ;r) = i^G([2rv];r) 



(9) 



and define Gcr(r) = iG([2no-]; r). Thus, for spin- 
polarized systems G(r) — J^a Ga{^), where 



G,(r) =aV 



w 



in which = \Vna-\'^ /8n„ and 



~ cx(conv) 



Note that ^ J2a 
spin-scaling relation 



cx(conv) 



(10) 



(11) 



because the 



applies only to energy densities. 



B. Construction by a coordinate transformation of 
the exact-exchange hole 

The exact-exchange energy density can be also con- 
verted to the gauge of a semilocal approximation by 



4 



transforming the exact-exchange hole. The conventional 
exact-exchange energy density can be written as 



„cx(conv) _ 



dr' 



where h^ai^^, r') is the exact-exchange hole 

l7a(r,r')P 



hxa{r, r') 



(12) 



(13) 



This hole is highly delocalized but can be made less 
so [H, m, [13] by an appropriate coordinate transforma- 
tion (r, r') — > (ri, r2) of the density matrix, such as [13] 



2 - cj -l + uj 

1 — LO LO 



(14) 



where < w < 1. This transformation does not affect the 
total exchange energy but yields a distinctly different 
exchange energy density 



(ri) 



"-^(ri) 



(15) 



where u = r2 — ri and /i"g.(ri, ri + u) is the transformed 
exact-exchange hole (HJ] 



/i^^(ri,ri +u) = K„{ri + [uj -l\\i,ri+uj\i) 

„ »cT(ri + l]u) 
"^(ri) 



(16) 



given in terms of the conventional exchange hole. Since 
the exchange hole associated with a semilocal functional 
is relatively local, a transformation of the exact-exchange 
hole by Eq. (fTB)) can make Cx^''"'' resemble its semilocal 
approximation more closely than e^'^°'^^'^ does. It should 



be noted that the transformed hole of Eq. p5)) does not 
obey the sum rule for the conventional hole at each ri, 
but preserves the correct normalization of the system- 
averaged exchange hole ^54?]. 

The extent of locality of the exchange hole depends 
on the value of parameter oj. The maximal localization 
is achieved at lo = 1/2 [U, [1^. We have numerically 
evaluated the transformed exact-exchange energy density 
exCT^'^'' (r) for various values of uj and found that uj = 0.92 
leads to the best fit of the exact-exchange energy density 
to the TPSS meta-GGA. 



III. COMPUTATIONAL METHODOLOGY 

In practice, it is much easier to construct the gauge cor- 
rection function G(r) than to perform numerical integra- 
tion over transformed coordinates of the exchange hole 
in Eq. (|15p . Therefore, we will adopt the former method 
for the purpose of constructing a hyper-GGA functional. 
In this section, we describe a general-purpose implemen- 
tation of the TPSS gauge term G{y) in finite basis sets. 



A. Evaluation of the exact-exchange energy 
density in the conventional gauge 

Analytic evaluation of the conventional exact-exchange 
energy density by Eqs. ^ and ^ is possible but im- 
practical because it requires evaluation and contraction 
of many one-electron integrals for each grid point r. In- 
stead, we employ a much more efficient approximate 
method of Delia Sala and Gorling [1^. Although this 
method is documented in the literature [2^, [1^ , we will 
supply its detailed derivation here because it serves as a 
stepping-stone for evaluating our function G{v). 

When a basis set {xti} is introduced, each Kohn-Sham 
orbital is taken as a linear combination of one-electron 
basis functions, <f)ia{v) = X^/i c^iX;i(i")- In terms of these 
basis functions, the density matrix of Eq. ([3]) is 



,(r,r')-^P;.X^(r)x:(r'), 



(17) 



where P^, = P^^ = E-' c^^^.)*- The conventional 
exact-exchange energy density of Eq. ^ can be written 
as 

Gx(conv)/ \ _ _f 



The single integral over r' in Eq. is not so easily 
evaluated for many different values of r, but introducing 
a second integration over r yields , which is evaluated 
analytically and simply in Gaussian basis sets. These 
facts motivate the following development. Using the S- 
function one can write 



Xr,(r) 



= dr" 



XrXr") 



(19) 



The (5-function can be approximated by an expansion in 
the same non-orthogonal basis as the orbitals, namely, 
S{r" — r) = ^^Xti{^)'^ti{^")j whose Fourier coefficients 
c^(r') can be determined as usual. This yields 



(20) 



where S^^ are matrix elements of the inverse of the basis 



set overlap matrix. Substitution of Eqs. p9|) and ([20] 
into Eq. (fTS]) gives 



gOx(co„v)(^) = i ^ ^ S;lKI^P^,X,{r)x:{r), (21) 
where 



(22) 



5 



are elements of the exchange matrix. Eq. ((2T|) can be 
rewritten as 



^cx(conv) 



i^) = \Y.^Ux,{v)xl{r), (23) 



where Q'J^^ are elements of the matrix Q"^ = S ^K°"P'^. 
Eq. ([23)) is analogous to the formula for the density 
n(r) EE 7<,(r,r) = E^^. ^/^i.Xp(r)x;^(r) except that, un- 
like P"", the matrix Q°" is generally not symmetric. The 
analogy can be made complete by replacing Q"^ with the 
symmetrized matrix 



1 



p<TJ^<Tg 1 ^ g Ij^ffp 



(24) 



The final formula for the conventional exact-exchange en- 
ergy density via the resolution of the identity is 

gCx(conv) (r) - 1 ^ Q^.Xm (r)x: (r) . (25) 

In practice, Cxct^'^""^^ is computed using the subroutines 
that evaluate n(r) by passing ^Q*^ in place of P'^. 



B. Evaluation of the exact-exchange energy density 
in the TPSS gauge 

The exact-exchange energy density in the TPSS gauge 
IS given by Eq. ©. The first term, ex^^™"''V), 

IS com- 
puted by Eq. (pS)) and the gauge term is evaluated as 
follows. Let us rewrite Eq. (fTU]) as 



G,(r) = a [V/,(r) • Ve,(r) + f^{v)V^e„{v)\ , (26) 
where 



/.(r) = 
Based on Eq. (fTTj). 



ric/e 



-2 



1 + 4c {n,/elf \ Ta 



Vcx(conv) ~ ^ 



(27) 



(28) 



.(29) 



Eqs. and involve the first and second derivatives 
of the exact-exchange energy density in the conventional 
gauge. These quantities are computed as 

^gCx(conv) ^ 1 ^ Q-^V[x^(r)x:(r)], (30) 
^2^ex(co„v) ^ 1 ^ Q^.V2[x4r)x:(r)], (31) 



0.00 


1 1 1 1 1 1 1 1 1 1 1 1 
Rco=2.132 bohr 


1.00 


111 




2.00 






3.00 


cc-pVTZ 






u-cc-pVTZ 




4.00 


u-6-311++G(3df,3pd) 




UGBS1P 

1 . 1 . 1 . 1 . 1 . 1 . 



-6 -4 -2 2 4 6 

z (bohr) 

FIG. 1: Conventional exact-exchange energy per electron in 
the CO molecule along the internuclear axis evaluated at the 
experimental geometry using the approximate resolution of 
the identity in four basis sets. The curves for the last three 
basis sets are close together everywhere. 
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FIG. 2: Exact-exchange energy per electron in the TPSS 
gauge in the NaCl molecule evaluated at the experimental 
geometry using the approximate resolution of the identity in 
four basis sets. The curves for the last three basis sets are 
almost indistinguishable. 



using the same subroutines that evaluate Vno-(r) and 
Verier (r) by passing ^Q"^ instead of P"^. Note that 
Eq. (|2ip clearly shows that matrix elements QJ^^, them- 
selves do not depend on r. 

The gradient V/o-(r) can be written as 



V/.(r) 
where 



V5,(r)+.g,(r)V 



5<T(r) 



(32) 



(33) 



el + 4cn2 ■ 

The quantity Vgo- (r) is evaluated using the chain rule as 
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usual and it involves only the first derivatives of n(r) and 
^cx(conv) ^ The second term can be written as 

The gradient involves derivatives of the type 

dx \ dx dx"^ dy dydx dz dzdx J 

(35) 

and similar expressions for n^^^ j dy and d\'S/na\'^ j dz. 
Finally, the gradient of the Kohn-Sham kinetic energy 
density Vto- has the components 

dx ^ \ dx dx'^ dy dydx dz dzdx J 

(36) 

and similarly for dr^/dy and dr^jdz. 

The quantities given by Eqs. ([35]) and ([36]) are not 
used in any of the common GGA and meta-GGA func- 
tionals and may not be immediately available in standard 
density functional codes. However, the first and second 
derivatives of the orbitals, from which Eqs. ([35| and ((36|) 
are built, are readily available. Thus, evaluation of the 
exact-exchange energy density and the gauge correction 
requires some modification of existing subroutines. We 
have implemented these formulas in a development ver- 
sion of the GAUSSIAN program [57| . 

C. Basis set effects 
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FIG. 3: Radial exchange energy densities of the H atom com- 
puted at the exact ground-state density: exact conventional 
[ex(conv)], exact in the TPSS gauge [ex(TPSS)], exact from 
a transformed exchange hole [ex(Lj), = 0.92], and TPSS. 
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We use the same nonorthogonal basis set {x^i} to ex- 
pand the Kohn-Sham orbitals and to approximate the 
(5- function by Eq. Since Eq. fSO]) in a finite ba- 

sis set is not exact, the conventional exact-exchange en- 
ergy density g"'^^'^""^) ^ylA its derivatives are only approx- 
imate when evaluated by Eqs. ([3D]), and (pij) . In 
fact, small and medium-size contracted basis sets may 
cause large errors in g^'^^'^""^) that are further magni- 

fied in er^^'''''^ via Ver^'°"'^ and V^e^^^^'^ Figs. [T] 
and[2]show that, for instance, the cc-pVTZ basis is insuf- 
ficiently flexible. On the other hand, the uncontracted cc- 
pVTZ basis set works almost as well as the near-complete 
UGBSIP basis [H, [5^. In general, uncontracted basis 
sets work much better in resolution of the identity tech- 
niques than the corresponding contracted bases. 

Furthermore, when cuspless Gaussian-type basis func- 
tions are used, |Vn(jp and IVtct^ exhibit spurious os- 
cillations in the vicinity of a nucleus. However, these 
artifacts are common to all semilocal density functional 
calculations employing Gaussian- type orbitals, are negli- 
gible energetically and may be ignored. 

In summary, we caution against using medium-size 
contracted basis sets like cc-pVTZ or 6-311-l-G* in 
Eqs. dUD, dSni), and When in doubt, it is always 



-0.8 I • ' • ' • ' • ' • 1 
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FIG. 4: Exchange energies per electron in the H atom com- 
puted at the exact ground-state density. For the explanation 
of the legend, refer to Fig. |3l 

safer to uncontract the basis set. In particular, we recom- 
mend the fully uncontracted 6-31H--|-G(3df,3pd) basis 
set, denoted as u-6-311++G(3df,3pd), which strikes per- 
fect balance between accuracy and computational cost. 



IV. RESULTS 

The fact that the TPSS meta-GGA was designed to re- 
cover many exact properties [1, jl^l of the exact-exchange 
functional does not guarantee that ej^^^ is close to 
^cx(conv)^ This is evident from Fig. [3] which shows radial 
plots of these energy densities in the H atom. The TPSS 
and conventional exact-exchange energy density are dif- 
ferent, even though they both integrate to the same exact 
value of —5/16 hartree Q. The exact-exchange energy 
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0.001 



FIG. 5: Radial exchange energy densities of the Ne atom com- 
puted at the converged Hartree-Fock orbitals in the UGBS 
basis set: exact conventional [ex(conv)], exact in the TPSS 
gauge [ex(TPSS)], and semilocal TPSS approximation. 




0.001 



0.01 0.1 
r (bohr) 



FIG. 6: Same as in Fig. [5]for the Kr atom. 



densities in the TPSS gauge ex'^^"'"^^^^ and ex'^^"-' are both 
much closer than g°'^('=°"^) to the semilocal ej^^^. 

Fig. m shows the exchange energy per electron, — 
Cx/n. vs. n for the H atom, comparing the exact con- 
ventional, exact in the TPSS gauge, and the hole- 
transformed {u! — 0.92) exact-exchange energies per elec- 
tron to the semilocal TPSS exchange approximation. Un- 
like Fig.[3l this figure shows what happens in the energeti- 
cally unimportant small-r and large-r regions. The trans- 
formed exact- exchange energy per electron Ex'^^"'' stands 
apart from the others in that it appears to have an in- 
verted cusp at the nucleus. 

Figs. [5] and [H] compare exchange energy in the Ne and 
Kr atoms evaluated in a post-self-consistent manner at 
the converged Hartree-Fock orbitals obtained using the 
near-complete UGBS basis set Overall, Figs. [3HB] 

suggest that, for the smaller atoms, the exact-exchange 






z (bohr) 



FIG. 7: Difference ej^^^— e™, where e^" is in the conventional 
and TPSS gauges, in a free N atom and along the internuciear 
axis of the N2 molecule at the experimental geometry. Panel 
b) shows only the right haif of the molecule with the N nu- 
cleus placed a.t z — 0. All quantities were computed at the 
converged TPSS orbitals using the UGBSIP basis set. 









b) 0.5 bohr off the axis 



RHeHe= 2-154 bohr 
RHeHe= 16bohr 



1-1 1 

z (bohr) 



FIG. 8: Difference ej^^^— e™''''^^^' in the He J molecule along 
the internuciear axis and along a parallel axis offset by 0.5 
bohr. Each panel shows only the region near the right nucleus 
which is always placed at 2 = 0. The static correlation in 
the stretched molecule He^'^^ • ■ ■ He^'^^ (dashed line) is more 
negative than at the equilibrium TPSS/cc-pVQZ geometry 
(solid line). All quantities were computed at the converged 
TPSS orbitals using the uncontracted cc-pVQZ basis set. 



energy density in the TPSS gauge is closer than g^'^^'^""^) 
to the TPSS exchange energy density. For larger atoms, 

T cx(TPSS) . 1 . ex(conv) . 

however, ex remains closer to Cx than to 

gTPSS the deep core region. 

Fig. [7] shows that the exact-exchange energy density 
(per electron) in the gauge of a semilocal approximation, 
^cx(TPSS)^ differs from jj^ a non-trivial way. The 

difference ej^^^ — gj^'''^"^^^^^ reveals subtle effects in the 
N2 molecule that are absent in a free N atom. Off-axis 
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Spherical iellium cluster 
r,=4, N=2 

ex(conv) 



ex(TPSS) 



TPSS 



6 8 10 
r (bohr) 



12 



14 



FIG. 9: Radial exchange energy densities for a spherical jel- 
lium cluster of A'^ = 2 electrons computed at the exchange- 
only OEP orbitals: exact conventional [ex(conv)] , exact in the 
TPSS gauge [ex(TPSS)], and semilocal (TPSS). 




Spherical ieiiium cluster 
rs=4, H=5h 

ex(conv) 



ex 

,ex(TPSS) 
TPSS 



10 15 20 
r (bohr) 



25 



30 



FIG. 10: Same as in Fig. [9]for a cluster of = 58 electrons. 



effects (not shown) are of course important. In contrast, 
the difference ej^^^ — g^'^^'^""^) jsj^ molecule is very 



ox(TPSS) 

Ex rep- 



similar to that in a free N atom. 

Fig. [5] shows that the difference ej^^® 
resenting the static correlation gets substantially more 
negative upon bond stretching, when the fragments show 
large fluctuations of electron number at the Hartree-Fock 
level, as it should. Note that the nuclei in the stretched 
HeJ molecule (RhcHo = 16 bohr) are essentially isolated, 
so the difference ej^^^ — g^'*^'"''^^^) jg almost perfectly 
symmetric about z = 0. 

We have also evaluated the exact-exchange energy den- 
sity in the TPSS gauge and the two conventional en- 
ergy densities for spherical jcUium clusters. A spheri- 
cal jellium cluster is a model system that has a uniform 
positive background charge and a spherically distributed 
electron density. The radius of the sphere is given by 



R = TsN^/^, where rg is the bulk density parameter and 
N is the number of electrons in the system. The volume 
of sphere is proportional to N and is given by the relation 
V = (47r/3)A^rj . The three exchange energy densities 
for the jellium clusters of = 4 for iV = 2 and 58 are 
shown in Figs. [51 and [TUl All three quantities, g^'^^'^""^)^ 
^cx(TPSS)^ and ej^^^, were evaluated at the orbitals and 
densities obtained from OEP calculations [H, [s^, [60j . 
As with atoms and molecules, Cx'^''"'"^^^^ is closer than 
to ej^^^ in these jellium clusters. 



cxfconv 

ex 



V. CONCLUSION 

As observed before [H, IH, |4^, the exchange energy 
density of a semilocal functional is reasonably close to the 
conventional exact-exchange energy density of Eq. ^ in 
compact systems like atoms or spherical jellium clusters. 
We confirm this here for the nonempirical TPSS meta- 
GGA. The relative differences are largest in regions of 
space where the density is dominated by a single orbital 
shape, making Jt close to 1, e.g., the H or He atoms 
and the two-electron jellium cluster. Particularly in these 
regions, the difference can be reduced by a gauge trans- 
formation of the conventional exact-exchange energy den- 
sity. 

We have found a simple, realistic, and not too highly 
parametrized form of the function G(r), given by Eq. ([7]), 
which via Eq. transforms the conventional difference 
of semilocal and exact-exchange energy densities appear- 
ing in Eq. (g]) to the gauge of the TPSS meta-GGA. 
This transformation solves the problem of nonunique- 
ness of the exact- exchange energy density arising in the 
context of modeling the static correlation by the differ- 
ence of semilocal and exact exchange energy densities. 

The transformed exact-exchange energy density Cx''*'"'"^^^'' 
does in fact contain more information about electron cor- 
relation than g™('^°'i'^)^ jji a, forthcoming article [4l[, we 
will present a construction of a hyper-GGA that relies on 
this gauge transformation to give highly accurate ther- 
mochemistry and reaction barriers. 

Finally, we have demonstrated that, as expected p^ . 
the difference between semilocal and exact-exchange en- 
ergy densities becomes more negative under bond stretch- 
ing in HeJ and related systems, where the separating 
fragments show large fluctuations of electron number at 
the independent-electron level. 
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